Magnetic quantum tunnelling in Fes with excited nuclei 
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We investigate the effect of dynamic nuclear spin fluctuation on quantum tunneling of the mag- 
netization (QTM) in the molecular magnet Fes by increasing the nuclei temperature using radio 
frequency (RF) pulses before the hysteresis loop measurements. The RF pulses do not change the 
electrons spin temperature. Independently we show that the nuclear spin-spin relaxation time T2 has 
strong temperature dependence. Nevertheless, we found no effect of the nuclear spin temperature 
on the tunneling probability. This suggests that in our experimental conditions only the hyperfine 
field strength is relevant for QTM. We demonstrate theoretically how this can occur. 

PACS numbers: 



in 

h 



I 

G 

O 

o 



> 

00 

q 
o 

On 
O 



X 



The importance of nuclei to quantum tunnelling of the 
magnetization (QTM) in Fes, was demonstrated exper- 
imentally by Wernsdorfer et al. [H, 0- They compared 
the tunnel splitting A of the standard Fes sample with 
a deuterated sample and a sample where 56 Fe was re- 
placed partially with 57 Fe. The enrichment with deu- 
terium causes a decrease of A, in accord with the de- 
creased hyperfine field (HF) Q. Similar conclusion was 
obtained with the 57 Fe enrichment. However, the ex- 
change of isotopes does not only vary the strength of the 
HF exerted on the molecule: it also changes the nuclear 
spin-spin relaxation rate T2. Both quantities might be 
important for the nuclear-assisted tunnelling process [H . 
Isotope substitution can not tell if only one or both quan- 
tities are relevant. Therefore, it is not yet established 
experimentally how exactly nuclei impact the tunnelling 
process. 

The experiment reported here aims at distinguishing 
between the contribution of the HF and T2. We present 
magnetization measurements of Fes during field sweep 
after transmitting radio frequency (RF) at the protons 
resonance. This transmission raises the protons temper- 
ature without changing the electrons temperature due 
to the enormous proton spin-lattice relaxation time T\ 
which is longer than 1000 sec at sub-Kelvin temperatures 
We also present proton spin-spin relaxation rate T2 
measurements and show that T2 varies with temperara- 
ture. When internal field fluctuations are slow, as in our 
case, T2 is a property internal to the nuclear spin system. 
Therefore, we argue that razing the nuclei temperature 
must lead to Ti variations with no modification to the 
hyperfine field. Our major finding is that T2 is not a 
relevant parameter for QTM. 

For this experiment, a Faraday force magnetometer 
shown in Fig. [1] was constructed inside the inner vac- 
uum chamber of a dilution refrigerator (DR) following 
the design of Sakakibara et al. [6(, with the addition of 
an RF coil. This magnetometer is suitable for measure- 
ments in high fields and at sub-Kelvin temperatures with 




FIG. 1: (Color online) Cross sectional view of the Faraday 
balance with: (1) movable plate of the capacitor, (2) screw 
for capacitor's fixed plate height adjustment, (3) sample, (4) 
PCTFE, (5) gold plated casing of the thermometer, (6) ther- 
mal link to the DR mixing chamber, (7) main coil, (8) gradient 
coils, (9) RF coil. 



no metallic parts near the sample. This is important for 
minimizing warming metallic parts with the RF. The DR 
is equipped with a main superconducting magnet that 
produces the field H, and two oppositely wound super- 
conducting magnets that produce a field gradient. 

The sample is mounted on the small load-sensing de- 
vice made of two parallel plates variable capacitor. The 
movable plate is suspended by two pairs of orthogo- 
nal crossed 0.2 mm diameter phosphor bronze wires at- 
tached to it with cpoxy. The static lower plate was 
mounted on an epoxy screw, for adjusting the initial 
capacity Co- When the sample is subjected to a spa- 
tially varying magnetic field B, it will experience a force 
F = M z (dB z / dz)z. This force is balanced by the wires. 
The displacement of the plate is proportional to F and 
can be detected as a capacitance C change. The total 
capacitance response is then given by 



Co 1 



C 



a-M z (dB z /dz) 



2 



where a is a constant that depends on the elastic prop- 
erties of the wires. 

The sample is grown by the method described in Rcf 
and is 20 mm 3 . It is oriented with its easy axis parallel 
to the magnetic field H . The sample is glued with GE- 
varnish to Poly-Chloro-TriFluoro-Ethylene (PCTFE), a 
fluorocarbon-based polymer, which has no hydrogen and 
is suitable for cryogenic applications. The bottom of 
the PCTFE is connected by a thermal link to the DR 
mixing chamber which produces the cooling, and to the 
movable plate. Approximately 2 cm above the sample, 
on the thermal link, there is a calibrated thermometer 
(RuC>2 R2200) in a gold plated casing. It is important 
to mention that the sample is in vacuum with no ex- 
change gas, and therefore its temperature T is not exactly 
the same as the temperature of the thermometer. How- 
ever, this is not a problem in our experiment since below 
400 mK the magnetization jumps of Fes are temperature- 
independent Q. 

In the magnetization experiments we apply a field of 
+1 T and wait until thermal equilibrium is reached. 
We then record the field value [Fig. [DJa)], capacitance 
[Fig. H[b)], and temperature [Fig. H^c)] as the field is 
swept from +1 T to -1 T at a rate dH/dt = 0.5 T/min. 
While we sweep the magnetic field from positive to neg- 
ative, we stop for several seconds at0.3 T (12.71 MHz) 
where we transmit the RF in the form of pulses as shown 
in Fig. [U(d). All attempts to deliver RF at a negative 
field resulted in immediate magnetization jumps, hence, 
the choice to transmit at a positive field. During the 
transmission, the temperature rises by 20 mK. When the 
field changes sign there is a larger temperature increase 
of 150 mK due to eddy currents. None of these tem- 
perature changes are enough to generate magnetization 
changes. To make the measurement with and without 
RF as similar as possible, we stopped at 0.3 T for several 
seconds even when we do not transmit RF. 

We first concentrate on the capacitance versus time, 
for a full sweep shown in Fig. [UJb). At times where the 
field is positive the capacitance is a smooth function of 
time (and field). This is because the spins are at their 
ground state for all positive fields and have nowhere to 
tunnel to. Once the field becomes negative, clear jumps 
in the capacitance are observed, indicating jumps in the 
magnetization that are taking place when tunnelling oc- 
curs between molecular spin states. The time it takes to 
sweep from the end of the RF transmission to the first 
jump is At = 60 sec. This time is much shorter than the 
nuclear Ti, as demonstrated in Fig. [2Jd). Therefore, the 
nuclei are expected to be excited when the Fes spins are 
tunnelling. We avoided RF transmission at fields higher 
than 0.3 T, and used high sweep rate in order to keep 
At short. Finally, Fig. [2]Jc) shows that magnetization 
jumps are accompanied by temperature spikes. These 
are discussed in a separate paper 

The results of measurements with and without the RF 
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FIG. 2: (Color online) The scheme of the measurements show- 
ing: the magnetic field sweept from positive to negative (a), 
the capacitance (b), the temperature (c), and the RF trans- 
mission (d). At is the time from the transmission to the first 
capacitance (magnetization) jump. Ti is the nuclear spin- 
lattice relaxation time. 

are summarized in Fig. [3J We focus on the first magneti- 
zation jump which is closest to the time of RF irradiation. 
The solid lines show sweeps with RF and the solid with 
symbols lines are sweeps without RF. We repeated these 
runs several times and found that within our experimen- 
tal resolution, and stability between individual sweeps, 
no effect of the RF can be detected. 

To appreciate this result we performed T2 measure- 
ments inside the mixing chamber of a DR and He cryo- 
stat using a more standard NMR setup and coil at fields 
of 0.76 T and 0.65 T and frequencies of 32 MHz and 
29 MHz respectivly. At these conditions the resonance 
field for most of the protons is not shifted from the free 
proton resonance, and the line width AH is on the or- 
der of 200 mT [lo| . The results of I/T2 are presented in 
the inset of Fig. [31 T2 varies from less than 10~ 4 sec at 
T = 3 K to 10~ 3 sec below T = 0.5 K. Between 4 K and 
150 K T2 is so short that no signal could be detected. 
Finally, due to the huge time scale difference between T2 
and T\ at all temperatures, it is reasonable to assume 
that Ti is determined by nuclear spin-spin coupling only. 

In the set up with both RF and magnetization shown 
in Fig. |T] it is difficult to detect the proton signal due to 
the poor filling factor in the Hclmholtz coil, the broad 
line width, and the extremely long T\. However, the 
GE- varnish gluing the sample has relatively narrow line 
and shorter T\. We therefore use the varnish signal at 



FIG. 3: (Color online) Capacitance measurements as a func- 
tion of field swept from positive to negative, with and without 
RF. The inset shows the temperature dependence of I/T2 in 
Fes measured near the free proton resonance condition. 

T = 140 mK to confirm the delivery of the RF radiation 
to the sample, to measure the strength of the RF field 
Hi, and to test our ability to saturate the nuclear tran- 
sitions. First, we measured the echo intensity as a func- 
tion of applied field at constant frequency of 12.71 MHz 
(0.3 T) using a n/2 — n pulse sequence. As shown in 
Fig. BJa), the full width at half maximum is only 4 ± 0.5 
mT. A similar line width was found for the varnish in our 
standard NMR spectrometer at 5 K. This ensures that we 
deliver the radiation to the center of the RF coil. Second, 
we determined the optimal pulse length. The echo inten- 
sity as a function of the pulse length f w / 2 is presented 
in Fig. 0Jb) . The maximum echo intensity was found at 
tir/2 =1-5 ± 0.5 fisec. From jHit^ = n/2 we calculated 
Hi to be 24 ± 4 mT. Finally, we determined T\, as pre- 
sented in Fig. rjjc) by saturating the proton transitions 
with a train of pulses, and then measuring the recovery 
of the signal at a time t using 7r/2 — n pulses. The pulse 
train equilibrates the up and down proton spins popula- 
tion. We found that the GE-varnish T\ is only 100 sec. 
More importantly, Fig. rjjc) demonstrates our ability to 
saturate the proton transitions. 

The above measurement allows us to estimate the vari- 
ation in the nuclear Ti at the time electronic spins are 
tunneling due to our RF irradiation. First, we exam- 
ine how many protons we excite. Since our H\ is smaller 
than the Feg line width AH, our direct pulses excite only 
Hi/ AH = 10% of the total number of protons. However 
during the transmission and after it, spin diffusion is tak- 
ing place spreading the nuclear temperature among all 
nuclei. The diffusion coefficient D is given by D = Wr 2 
where W is flip-flop rate of neighboring nuclei, and r 
is the distance between them [11]. For dipolar coupling 
W > (~f 2 h/r 3 ) 2 /(jAH) where j 2 h/r 3 is the strength of 
the dipolar interaction, and 1/jAH is a lower limit on 
the density of states [llj. The time it takes for the heat 
to spred among all nuclei in a unite cell of volume V is 
V 2 / 3, /D, which is less than 10 sec. Therefore, all nuclei 
should be warm before the first tunnelling event is taking 
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FIG. 4: (Color online) Echo intensity at 140 mK from the GE- 
varnish gluing the sample as a function of field (in frequency 
units) (a), pulse length (b), and time after saturation (c). The 
solid lines are guides to the eye. 

place. 

Second, we evaluate by how much the irradiate nu- 
clei cool during the time between transmission and tun- 
nelling. For this we employ the equation 

i-(l-exp(-At/Tx)) = J- 

-*- e n 

where T e is the electron's temperature (140 mK) and 
T n is the nuclei temperature. Immediately after the RF 
pulses (At = 0) T n = 00. As At grows T„ decreases until 
at At = 00 it reaches the electrons temperature again. 
This equation suggests that the nuclei temperature at the 
time of the tunnelling is well above 3 K where T2 increases 
by a factor 14 from its value at 140 mK. Therefore, we 
conclude that changing T2 by an order of magnitude has 
no effect on the tunnelling probability, for our sweep rate. 

To demonstrate that it is conceivable to have an iso- 
tope effect, yet no dependence on T2, we analyze an effec- 
tive model for the dynamics of the system in the vicinity 
of a resonant transition between molecular spin levels m 
and to'. This is essentially the Landau-Zencr (LZ) prob- 
lem with the addition of a transverse magnetic noise. 
The effective Hamiltonian, describing a spin 1/2 with a 
resonance tunnel splitting A, a time-dependent magnetic 
field in the z direction, and a fluctuating magnetic field 
in the x direction, is given by 

H = atS z + AS X + B%(t)S x ■ (1) 

Here S x = \a x and S z = \o z , where a x and a z are 
the Pauli matrices, and a = 2gfiB (wi — Tn') dH/dt. We 
assume that the stochastic field B%(t) has a correlation 
function 

{Bl(t)B e x (t')) = (Bf) exp(- |f - f'| /r c ) (2) 
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where r c is the correlation time. B e x is related to the 
hypcrfine field (see below), and r c stands for T^. We 
consider only a transverse fluctuating field since for the 
— 10 to 9 transition, the measured tunnel splitting A = 
3 x 10~ 7 K [l2j|, and our sweep rate, the sudden limit is 
obeyed, namely, A/^/Ka < 1. In this case it is well es- 
tablished that a stochastic field coupled to the z direction 
of the spin has no effect on the LZ tunnelling probability 

We next write the wave function as &(t) = C_ (t) \ — ) + 
C-f(i) |+), where |±) denote eigenstates of a z . Defining 
C±(t) = cxp(±iat 2 /Ah)C±(t) and introducing a dimcn- 
sionless time variable y = t^Ja/fi, the Schrodinger equa- 
tion can be expressed in the integral form 



C T (oo) 



>±iy 2 /2 



: f°° (A + B e x (y)) 

J — OO 



(3) 



C±(y)dy. 



Assuming the initial conditions C+(— oo) = 1 and 
C_(— oo) = 0, the tunnelling probability is given by 

P = 1 — ^|C_(oo)| 2 ^) where () stands for an average 

of stochastic field realizations. In the sudden limit, C+ 
does not change much. If, in addition, B x -C V ah then 
C+ under the integral can be replaced by 1 to first order 
in A and B. x . This yields 



P= 1- 



1 

4ak 



dxdy (A 2 + (Bf) e^~y\y<y 2 ' x2 )/ 2 



(4) 

where we have used Eq. [2] with v = l/r c y/ a/h. This 
yields 



P=l-n(A 2 + (Bf))/2ha, 



(5) 



in which there is no dependence on the parameter v. The 
transition probability is therefore dependent on the HF 
strength, but not on its correlation time r c . 

1/2 

To make this derivation attractive {B x 2 ) must be 
of the order of the measured tunnel splitting. When con- 
verting the Fes problem to the two level LZ problem B x 
is scaled down from the field B x the nuclei produce, since 
B x has a matrix element between the m and m! states 
only in the |m — m'|'th order of perturbation theory. As a 
consequence, Garanin and Chudnovsky [Til ] showed that 



Bt 



2D 



1 (S + m')\(S - m)\ 



(m' - m - l)! 2 V (S - m')\(S + m)\ 



B x 
2D 



(6) 

where D = 0.27 K is the Fes single ion anisotropy coeffi- 
cient. Protons produce a field of the order of 1 — 10 mT 
inside a solid, corresponding to B x of 0.01 — 0.001 K, 



which is not small. However, in our case m! 
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there has to be a shortcut in the tunneling process that 
will make m' — m smaller and of the order of 6 — 7. It 
is reasonable that such a shortcut exists since we, and 
other researchers 12j|, see only four magnetization jumps 



and not 10. 

To summarize, we exploit the strong temperature de- 
pendence of the nuclear spin-spin relaxation time T2 
around 1 K in order to test the effect of nuclear fluctua- 
tions on quantum tunneling of the magnetization. Since 
in our case T2 is a property internal to the nuclear spin 
system, we change it by warming only this system with 
radio frequency transmitted at the protons resonance. 
We then measure the size of the magnetization jumps 
due to tunneling. During the magnetization measure- 
ments the nuclei stay warm due to the enormously long 
spin-lattice relaxation time T\. We found no effect of the 
nuclear spin temperature on the magnetization jump and 
argue that T2 is a parameter irrelevant to the tunneling 
probability in our experimental conditions. We present a 
calculation demonstrating that nuclear spins can, indeed, 
affect the tunneling via their hypcrfine field strength only. 
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this research. 
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